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the characteristic mantle turnover time is R/u = 58 Myr. Thus we find that x 
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Geoid, Topography, and Convection-Driven Crustal Deformation on Venus; 

Mark Simons, Bradford H. Hager, and Sean C. Solomon, Department of Earth, Atmospheric, 
and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139. 

Introduction. High-resolution Magellan images and altimetry of Venus reveal a wide 
range of styles and scales of surface deformation [1] that cannot readily be explained within 
the classical terrestrial plate tectonic paradigm. The high correlation of long-wavelength 
topography and gravity and the large apparent depths of compensation suggest that Venus 
lacks an upper-mantle low- viscosity zone [2-5]. A key difference between Earth and Venus 
may be the degree of coupling between the convecting mantle and the overlying lithosphere. 
Mantle flow should then have recognizable signatures in the relationships between surface 
topography, crustal deformation, and the observed gravity field [6,7]. 

Model. We explore the effects of this coupling by means of a finite element modelling 
technique. The crust and mantle in these models are treated as viscous fluids. We solve 
both the equations of motion and the heat equation at every time step using a modified 
version of the 2-D Cartesian finite-element program ConMan [8]. A passive marker chain 
tracks the crust-mantle interface and permits variation in the crustal buoyancy as well as 
specific crustal and mantle rheologies. These rheologies depend on composition, temperature 
and stress. In addition to the flow field, the stress field in the lithosphere, the surface 
topography, and the resulting geoid are readily calculated. The models presented here use an 
irregular finite-element mesh that is 50 elements high and 160 elements wide. Our maximum 
resolution is in the 40-km-thick top layer, where each element is 2 km high and 5 km wide. In 
all, the mesh is 800 km in the horizontal dimension and 400 km in the vertical dimension. We 
impose free-slip boundary conditions on the top and side walls, with no flow through these 
walls. Flow at the bottom boundary is constrained to be vertical with no horizontal flow 
permitted. This last boundary condition gives a virtual 800 km by 800 km box. The surface 
topography is calculated from the vertical stresses on the top wall of the box. Top and bottom 
temperatures are fixed at 500°C and 1250°C, respectively. Initially, we impose a linear 
temperature gradient across the lithosphere and set the rest of the mantle to be isothermal. 
We investigate two classes of models. Flow in the first class is initiated with a sinusoi dal 
temperature perturbation throughout the box. This results in mantle flow that is dominated 
by concentrated down welling (Figure lb). In the second class of models, flow is driven by a 
hot patch at the side of the box, thereby driving flow by concentrated up welling (Figure la). 



Horizontal Distance, km Horizontal Distance, km 


Figure 1. Temperature contours at 130 My for the (a) upwelling and (b) down welling models. 
Contour interval is 100 °C; the top of the box is at 500 °C. 
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Results. In all our models, convection produces horizontal compressional stresses in 
lithosphere above downwelling mantle and extensional stresses in lithosphere above 
upwelling mantle. As the convective vigor increases so does the magnitude of the stress. In 
models with constant-viscosity mantle overlain by a constant-viscosity crust, stress in the 
crust reaches values in excess of 100 MPa in less than 100 My. We find that the rate of 
increase in compressive stress decreases with increasing crustal viscosity. This is because 
the stronger the crust, the more the development of the convective instability in the mantle 
driving the deformation is impeded. We also find that the magnitude of the peak compressive 
stress achieved above the downwelling increases with higher viscosities and/or with thinner 
initial crustal layers; the stronger the crustal lid, the more are tractions from mantle 
convection supported in the crust. Since force balance on the crust requires that shear 
traction integrated along the base be balanced by normal tractions integrated through its 
thickness, the thinner the crust, the larger the horizontal stresses. 

Both analytical models [9,10] and our numerical models of convection-induced crustal flow 
indicate that the amplitude and sign of the topography are highly time- and rheology- 
dependent. In general, possible responses of the crust to mantle flow can be divided into 
three categories. The first involves little, if any, crustal flow, and topography results mainly 
from the transmission of normal tractions induced by density contrasts within the mantle. 

The second possible regime involves substantial crustal flow, with geologically rapid 
thickening over convective downwelling and thinning over convective upwelling. In this 
regime the effects of crustal thickness variations dominate the topography. A third possible 
regime lies between the first two, with "in phase" deformation on short time scales and 
crustal flow on longer time scales. A strong mantle lithosphere tends to shield the crust from 
convective shear tractions, and topography results mainly from the transmission of convective 
normal tractions. A relatively weak lower crust facilitates crustal deformation, and the 
isostatic effects of crustal thickness variations dominate the topography. 

Consideration of geoid to topography ratios (GTRs) can restrict which regime of crustal 
response is appropriate for Venus. The distibution of estimated GTRs for several highland 
regions on Venus is bimodal with two clusters around 10 and 25 m/km [5]. The positive 
correlation of long-wavelength gravity and topography implies that there are no major regions 
that have negative GTRs. To keep the GTRs positive, the geoid must follow the surface 
topogaphy in sign at all times. Although the regime of negligible crustal flow can correctly 
predict both the sign and magnitude of the GTRs, it does not allow for the crustal deformation 
(i.e., flow) inferred from observations of tectonic features on Venus. 

In contrast to the regime of negligible flow, that of time-dependent crustal flow generates 
topography that changes sign. Over a mantle downwelling, the topography is negative in the 
early stages of deformation and positive in the later stages of deformation; the converse holds 
over a mantle upwelling. During the transitional period the topography goes through zero, the 
geoid does not go through zero, and the GTR is unbounded; this singularity is not observed 
on Venus. An example of the evolution of the GTR for a model in the time-dependent regime 
is shown in Figure 2. Note that before 90 My both topography and geoid are negative, but at 
90 My crustal thickening effects begin to dominate and the topography switches sign, causing 
the GTRs to first become unbounded and then negative. 

In the absence of a mechanism by which the sign of the geoid anomaly mimics that of the 
topography over a given upwelling or downwelling, only the regime of rapid crustal flow is 
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plausible. In the case of mantle downwelling this would also require that the lower mantle be 
more viscous than the upper mantle in order to produce the required positive geoid anomalies. 
This has already been shown to be true for the Earth, where the observed geoid highs over 
regions of mantle upwelling and regions of mantle downwelling are best explained by the 
presence of a strong lower mantle [11,12]. The large positive GTRs and the presence of 
large shield volcanoes in certain highland regions on Venus, such as Beta Regio and Eistla 
Regio, are best explained as areas of mantle, upwelling [5,13,14]. The regime of rapid crustal 
flow predicts crustal thinning over the upwelling. However, the extensive partial melt and 
ensuing volcanism expected over such regions of mantle may outweigh the effects of crustal 
thinning on the surface topography and thus also yield postive GTRs [15]. 



Figure 2. Evolution of GTR with time. At each time step, the GTR shown is the slope of the 
best fit line to the collection of geoid and topography values upward continued to 250 km 
above the surface nodes in the model. The result is for a downwelling model, with a crust of 
constant viscosity 50 times that of the mantle. In this model the sign of the topography is 
time-dependent, but that of the geoid is not. The singularity at about 90 My corresponds to 
the time when topography crosses the zero datum. Note that the GTRs for this model are 
neither always positive nor always bounded. 
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Introduction. High-resolution Magellan images and altimetry of Venus reveal a wide range of 
styles and scales of surface deformation [1] that cannot readily be explained within the classical 
terrestrial plate tectonic paradigm. The high correlation of long-wavelength topography and gravity 
and the large apparent depths of compensation suggest that Venus lads an upper-mantle low- v 
viscosity zone [2-5]. A key difference between Earth and Venus may be the degree of ooupling 
between the convecting mantle and the overlying lithosphere. Mantle flow should then have 
recognizable signatures in the relationships between the observed surface topography, crustal 
deformation, and the gravity field [6,7]. Therefore, comparison of model results with 
observational data qan help to constrain such parameters as crustal and thermal boundary layer 
thicknesses as well as the character of mantle flow below different Venusian features. We explore 
in this paper the effects of this coupling by means of a finite element modelling technique. 

Model The oust and mantle in the models are treated as viscous fluids'. We solve both the 
equations of motion and the heat equation at each time step using a modified version of the two- 
dimensional, Cartesian finite-element program ConMan [81. Our formulation assigns material 
properties to each element and explicitly tracks the free surface topography and the crust-mantle 
interface bjr letting the grid deform with time in a semi-Lagrangian fashion. This procedure 
permits variation in the crustal buoyancy as well as specific crustal and mantle rheologies that can 
depend on composition, temperature and stress. In addition to the flow field, the stress field in the 
lithosphere, the surface topography, and the resulting geoid are readily calculated. The model 
domain is a square box. We impose frec-slip boundary conditions on the bottom and side walls, 
with no flow through these walls. The top of the box is a true free surface, so there is no need to 
derive surface topography from vertical stresses on the top of the box. We investigate two classes 
of models. In the first class flow is dominated by concentrated upwelling (Figure la), and in the 
second class flow is dominated by concentrated downwelling (Figure lb). We vary the initial 
crustal and thermal boundaiy layer thicknesses as well as the effect of crustal and mantle viscosities 
that are either constant, temperature-dependent, or fully non-Newtonian. 

Results. In all of our models, convection produces horizontal compressional stresses in 
lithosphere above downwelling mantle and extensional stresses in lithosphere above upwelling 
mantle. As the convective vigor increases so does the magnitude of the stress. In models with 
constant-viscosity mantle overlain by a constant-viscosity oust, stress in the crust reaches values 
in excess of 100 MPa in less than 100 My. We find that the rate of increase in compressive stress 
decreases with increasing crustal viscosity. This is because the stronger the crust, the more the 
development of the convective instability in the mantle driving the deformation is impeded. We 
also find that the magnitude of the peak compressive stress achieved above the downwelling 
increases with higher viscosities and/or with thinner initial crustal layers; the stronger the crustal 
lid, the more are tractions from mantle convection supported in the crust. Since force balance on 
the crust requires that shear tractions integrated along the base be balanced by normal tractions 
integrated through its thickness, the thinner the crust, the larger the horizontal stresses. 

Both analytical models [9,10] and our numerical models of convection-induced crustal flow 
indicate that die amplitude and sign of the topography are time- and rheology-dependent. In 
general possible responses of the crust to mantle flow can be divided into three categories. The 
first involves little, if any, crustal flow, and topography results mainly from the transmission of 
normal tractions induced by density contrasts within the mantle. The second possible regime 
involves substantial crustal flow, with geologically rapid thickening over convective downwelling 
and thinning over convective upwelling. In this regime the effects of crustal thickness variations 
dominate the topography. A third possible regime lies between the first two, with "in phase" 
deformation on short time scales and crustal flow on longer time scales. A strong mantle 
lithosphere tends to shield the crust from convective shear tractions, and topography results mainly 
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